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Abstract 

By a closer inspection of the massive Schwinger model within mass perturbation theory 
we find that, in addition to the n-boson bound states, a further type of hybrid bound 
states has to be included into the model. Further we explicitly compute the decay widths 
of the three-boson bound state and of the lightest hybrid bound state. 
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1 Introduction 



The massive Schwinger model is two-dimensional QED with one massive fermion. In 
this model there are instanton-like gauge field configurations present, and, therefore, a 9 
vacuum has to be introduced as a new, physical vacuum ([|l], §])• Further, confinement is 
realized in this model in the sense that there are no fermions in the physical spectrum 
(P, H). The fermions form charge neutral bosons, and only the latter ones exist as 
physical particles. The fundamental particle of the theory is a massive, interacting boson 
with mass ^ = Mi (Schwinger boson). In addition, there exist ri-boson bound states. 
The two-boson bound state is stable (mass M2), whereas the higher bound states may 
decay into Mi and M2 particles ( [0 ) . AH these features have been discussed in |^ within 
mass perturbation theory (0, - which uses the exacly soluble massless Schwinger 
model ([0], - |]T3[) as a starting point. 

Here we will find that another type of unstable bound states has to be included into 
the theory, namely hybrid bound states composed of Mi and M2 particles. In addition, 
we will compute the decay widths of the M3 bound state and of the lightest hybrid bound 
state (which consists of one Mi and one M2 and has mass Mi_i). 



2 Bound states 

For later convenience we define the functions 

E±{x) = e^^"^''^") - 1 (1) 

and their Fourier transforms E±{p), where D^{x) is the massive scalar propagator. As 
was discussed in all the n-boson bound state masses M„ may be inferred from the 
two-point function (P = ^^75^, S = S± = ^^^l ± 75)^) 

n(a;) := 6{x) + g{P{x)P{0)) (2) 

(or g{S{x)S{0)) for even bound states) where g = mS -|- o(m^) is the coupling constant 
of the mass perturbation theory for vanishing vacuum angle 6 = (the general 9 case we 
discuss in a moment); E is the fermion condensate of the massless model. Il{x) is related 
to the bosonic n-point functions of the theory via the Dyson-Schwinger equations P]. In 
momentum space Il{p) may be resummed, 

flip) = i (3) 

where n.f. means non-factorizable and denotes all Feynman graphs that may not be 
factorized in momentum space. In lowest order (PP)^ ^ is 

{PP)nAp) = liEM-E-{p)) (4) 
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(and with a + for {SS)^f ). Expanding the exponential one finds 1 — gJ2'^=i ^;^^Jlip) 
in the denominator of (3) (more precisely, the odd powers for (PP)^ ^ , the even powers 
for {SS)^^ ). At p"^ = {nnY, D^{p) is singular, therefore there are mass poles = 
slightly below the n-boson thresholds. Further D^{p) have imaginary parts at p^ = for 
m < n, therefore decays into mMi are possible (more precisely, for the parity conserving 
case ^ = 0, only odd odd or even — even decays are possible). 

Up to now we did not mention the M2 particle, although decays into some M2 are 
perfectly possible. So where is it? The boson bound states are found by a resummation, 
so a further resummation^is a reasonable idea. Let us look at the Mi + M2 final state for 
definiteness. Within g{PP)^^ {p) we may find the following term 

Hip) := / j^g{m,AQ)n{q)g{PP)^, iq)A7rD,{p - q). (5) 

It is simply a loop where we have selected one boson to run along the first line, whereas 
the g{PP)^^ (q) and n(g) run along the other line. H{p) is a loop, therefore it is non- 
factorizable. The additional {g{PP)^^{q)Y factor is necessary, because n(g) starts at 
zeroth order (n(g) = 1 + g{PP)^^ (g) + . . .), and without this factor we would include 
some diagrams into H{p) that were already used for the n-boson bound-state formation 
(double counting). 

The claim is that H{p) has a threshold singularity precisely at p^ = (Mi + M2)^, and 
therefore an imaginary part for jo^ > (Mi + M2)^. But this is easy to see. At = M| 
n(g) has the M2 one-particle singularity and g{PP)^^ (M2) = 1. Therefore, near p^ = 
(Ml + M2)^, H{p) is just the Mi, M2- two-boson loop (up to a normalization constant). 

Observe that this line of reasoning is not true for higher bound states, p"^ ~ (M„ + 
Mi)^, n > 2. n(M„) contains imaginary parts and is not singular for n > 2 (because the 
M„ are unstable), and therefore H{p) has no thresholds at higher p'^. 

A further consequence is that H{p) gives rise to a further mass pole slightly below 
p2 = (Mi + M2)2 in (3). 

These considerations may be generalized, and we find niMi + n2M2 particle- 
production thresholds at = (riiMi + n2M2)^ and (unstable) niMi + n2A/f2-bound 
states slightly below. 

After all, this is not so surprizing. The M2 are stable particles and interacting via an 
attractive force. In two dimensions this must give rise to a bound state formation. (Similar 
conclusions may be drawn from unitarity when M2-scattering is considered, [0.) 

Before starting the actual computations, we should generalize to arbitrary 6^0. 
There the coupling constant is complex, g — > ge.gl, and, because of parity violation, the 
Feynman rules acquire a matrix structure (the propagators are 2x2 matrices, the vertices 
tensors, etc.). The exact propagator may be inverted, analogously to (3), and leads to 
(see i) 

^ (6) 

!-«-«* + aa* - (3(3* 

a{p) = ge{S+S+)^jXp) > P{p) = 9e{S+S_)^^i\p) (7) 
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and J^ij = +, — ) gives the {SiSj) component of the propagator. For our considera- 
tions only the denominator in (6) is important. In leading order 

a{p)=geE+{p) , P{p) = geE^p) , 9e = ^^e^' (8) 
and the denominator reads 

1 - mJ: cos 9 E^p) + '^f^{El{p) - El{p)). (9) 

Inserting the n-boson functions {dn{p) '■= (p) ) results in 

1 - mJ:cos9{di + d2 + ...) + m^Ts^(di{d2 + c/4 + . . .) + 4(^2 + ^4 + ...) + .. .) (10) 

Now suppose we are e.g. at the M3 bound state mass. Then the real part of (10) vanishes 
by definition and rnLcos6d^{M^) = 1 + o(m), and we get 

-imH cos 9\md2{M:i)+im^T?d^{M^)\md2{M^) = -zmS(cos^ ^)lmd2{M3). (11) 

cos 

This computation may be generalized easily, and we find that each parity allowed decay 
acquires a cos 6*, whereas a parity forbidden decay acquires a (cos 6* — ^^) factor. 

To include the decays into M2 we have to perform a further resummation analogous to 
above, however, the resummed contributions enter into the functions a, /5 in a way that 
is perfectly consistent with our parity considerations (a riiMi + n2M2-state has parity 

P = (-1)"!). 



3 Bound state masses 

We are now prepared for explicit computations, but before computing decay widths we 
need the masses and residues of the propagator at the various mass poles. The masses 
Ml, M2, M3 have already been computed (0; there is, however, a numerical error in the 
M2 mass formula in |^), 

Ml = = ^0 + Ai + o(m^) , Ai = 47rmS cos 9 (12) 
M| = 4/i2-A2 , A2 = 13 

Ml = ^l?~^-i , A3 ^ 6.993/i2exp(-0.263— f -) (14) 

mE cos 9 

and the three-boson binding energy is smaller than polynomial in the coupling constant 
m (or g). 

In leading order the n-th mass pole is the zero of the function 

Up'') = I- mT. cos 9dn{p'), (15) 
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2 A/f2\ 



therefore the residue may be inferred from the first Taylor coefficient around {p — M, 

/„(p2)^c„b2-M„2). (16) 

The Cn may be inferred from the computation of the mass pofes (0) and are related to 
the binding energies. Explicitly they read 

= 1 = !_ (17) 

47rmScos^ Ai ^ ' 

fi^ 1 

C2 



87r4(mScos^)2 2A2 ' ' 

mS cos 6 

^ 0.263//2A3 ^^^^ 

The mass Mi.i is the solution of 1 = {g0 + gg)H{p), which looks difficult to solve. However, 
there is an approximation. At threshold Il{q) equals the M2 propagator, so this may 
be a reasonable approximation provided that the binding energy is sufficiently small, 
Ai^i = (Ml + M2)^ — < A2. In this approximation we have for Mi^i 



1 = mS cos 9 



(Pq Syr^mS cos 6 An 
(27r)2 /i2(g2 _Ml){p-qy- Ml 
327r5m2s2 cos^ 9 



2^ 



7r + 



2Tiij?w{p^,MlMl 

2p^ \ 
^'"^^^ w{p\ Ml Ml) - {p' + Ml - M|)(p2 - + Ml) ) 

tt;(x, y, z) := (— — — + 2xy + 2x2; + 21/2)2 (21) 

where we inserted the residues that may be derived from the Taylor coefficients Ci , C2 
(17,18) (ReSj = — ^ — h)- The solution is 

\ ^ ' \ ' amS COS f) ' 

..2 ...2 A . 327riO(mScos^)^ 

M^ = Mi + M2f - Ai,i , Ai.i = ^— ^ (22) 

which shows that our approximation is justified for sufficiently small m. 

Mi^i was computed in a way analogous to M2 (see 0), therefore it leads to an anal- 
ogous Taylor coefficient 

1 fj,^ 
"^''^ " 2A^ " 647rio(mScos^)4' ^^^^ 
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4 Decay width computation 



The decay widths may be inferred in a simple way from the imaginary parts of the 
propagator. Generally 

r^i \ const. 

and r is the decay width. In our case the poles have their Taylor coefficients, 
~ const. const'. 

"^^^ ^ c.(p2-Mf)-zIm(---) ^ p2 _ Ml - i'-^ ^ ' 

and therefore 

- (26) 

Before performing the explicit computations let us add a short remark. The q are related 
to the binding energies, Cj ~ 7^-. Therefore, all the decay widths are restricted by the 
binding energies, Fj ~ Aj. But this is a very reasonable result. The denominator of 
the propagator (25) has zero real part at Mf and infinite real part at the real particle 
production threshold. Suppose Tl{p) contributes to a scattering process (to be discussed 
in detail in a further publication, []T3|). It will give rise to a local maximum (resonance) at 
= M^, and to a local minimum at the production threshold p"^ = Mf + Aj. Therefore 
the resonance width (decay width) must be bounded by Aj. 

Now let us perform the explicit calculations. At Ml ^ the propagator is 

^^^"^ ci,i(p2 - Mf 1) - zmS(cos6' - ■^)\md2{p) ^'^'^^ 



Im d2{p) 



2w{p'^,MlMl) 

w{x,y, z) = [x +y +z — 2xy — 2xz — 2yz)2 (28) 
leading to the decay width (Mi = fi) 

2V2(mScos^)^ 1 ^, ^^^^^ ,mcos6',., 1 
9v5/i^ cos^ n cos^ 

(S = ^ = 0.283/i) for the decay Mi 1 — > 2Mi. This decay is parity forbidden, and 
therefore Mi 1 is stable for 6 = 0. 

For the M3 decay there exist two channels, M3 M2 + Mi, M3 ^ 2Mi, 

"^^^ ^ csip' - Ml) - zmS(cosg - ^) ,(,.f4,^.) - ^(^^ cos g)^ ^.^^;^;^, ' 
leading to the partial decay widths 

rM3^2M. = 0.263^^(^ - 1) ^ 3.608//(^ - 1) exp(-0.929^^) (31) 
9v5^ cos^ 6* cos"^ t/ mcosO 



' — ^ ~ 43.9/iexp -0.929 — — 
3v3^ mcosO 



'm,^m,+m, = 0.263-^ ~ 43.9/i exp(-0.929^— ) (32) 



and to the ratio 



_J. 1 

M3^2Mi _ cos^e ^ (^33^ 



rM3^M2+Mi vTstt 

The latter is independent of the approximations that were used in the computation of 
M3 and C3. Observe that Tm3^M2+Mi is larger than rAf3^2Afi, although Mi + M2 ~ M3. 
This is so because the phase space "volume" does not rise with increasing momentum in 
d=l + l. 

Remark: there seems to be a cheating concerning the sign of Tm3^2Mi (see (30), (31)). 
This is a remnant of the Euclidean conventions that are implizit in our computations 
(see e.g. 0). There the conventions are such that 9 is imaginary and, consequently, 
cos^ — > 0. In a really Minkowskian computation, roughly speaking, the roles of 

and E_ in (6) are exchanged, leading to a relative sign between odd and even states. 
The final results (29), (31) and (32) are expressed for Minkowski space and for real 9 
(cos^ — 1 > 0), which explains the sign. 



5 Summary 

By a closer inspection of the massive Schwinger model we have found that its spectrum 
is richer than expected earlier. In addition to the n-boson bound states there exist hybrid 
bound states that are composed of fundamental bosons and stable two-boson bound 
states. A posteriori their existence is not too surprizing and may be traced back to the 
fact that particles that attract each other form at least one bound state in ci = 2; or 
it may be understood by some unitarity arguments. For the special case of vanishing 
vacuum angle, 9 = 0, the lowest of these hybrid bound states is even stable and must be 
added to the physical particles of the theory. 

Further we computed the decay widths of some unstable bound states and found that 
our results are consistent with an interpretation of the bound states as resonances. Even 
more insight into these features would be possible by a discussion of scattering, which 
will be done in a forthcoming publication (|T^]). 

Of course, it would be interesting to compare our results to other approaches, like e.g. 
lattice calculations. 
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